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ABSTRACT 

We present exact results for a lattice model of cluster growth in ID. The growth mech- 
anism involves interface hopping and pairwise annihilation supplemented by spontaneous 
creation of the stable-phase, +1, regions by overturning the unstable- phase, —1, spins with 
probability p. For cluster coarsening at phase coexistence, p = 0, the conventional structure- 
factor scaling applies. In this limit our model falls in the class of diffusion-limited reactions 
A+A— >inert. The +1 cluster size grows diffusively, ~ \/t, and the two-point correlation 
function obeys scaling. However, for p > 0, i.e., for the dynamics of formation of stable 
phase from unstable phase, we find that structure-factor scaling breaks down; the length 
scale associated with the size of the growing +1 clusters reflects only the short-distance 
properties of the two-point correlations. 



1. INTRODUCTION 

Lattice cellular automaton-type models with local tendency for ordering, termed voter 
models, can be used to study phase segregation and cluster coarsening reminiscent of spin- 
odal decomposition/ 1 ' 2 ^ at least in low dimensions. Both the cluster-size^ 3 -* and structure- 
factor*- 4 ) scaling at phase separation have been subjects of numerous investigations. How- 
ever, most of the available results for realistic 2D and 3D dynamical models are numerical. 
We distinguish between the two "scaling" terms as follows. By cluster-size scaling we mean 
scaling properties of the cluster size distribution. The term structure-factor scaling is re- 
served for the scaling properties of the two-point order parameter correlation function. The 
latter is accessible to scattering experiments. 
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The symmetric voter- type models are related also to the diffusion-limited chemical reac- 
tions involving particle annihilation, A+A^inert.^ 5 ) There are several exact results available 
mostly in ID which essentially translate to various average and asymptotic properties of 
the cluster size distribution in the phase-separation nomenclature. ( 5_14 ) Recent works also 
yielded exact results for the two-point correlations. ( 13 > 14 ) Results for the chemical reaction 
models in D > 1 are more limited: see Ref. 10 and literature cited therein. Furthermore, 
the relation of chemical reaction systems to voter models is less straightforward.* 1 ' 15 ~ 17 ) 

The purpose of the present work is to introduce a lattice model that incorporates voter- 
type cluster coarsening by interface diffusion in ID, as well as the process of spontaneous 
formation of stable-phase regions from those of the unstable phase. We derive exact results 
for the two-point correlations. Our main finding is that structure-factor scaling ideas cannot 
be extended from cluster coarsening (of both phases) at coexistence to stable-phase cluster 
growth off coexistence . While the general- D formulation is outlined, the present study is 
focused on the ID case. 

The model is defined in Section 2. The generating function solution of the discrete-time 
and discrete-space dynamics is presented in Section 3. Section 4 is devoted to the discus- 
sion of some special limits including the symmetric case. Our results are consistent with 
previous studies; the general framework of our formulation is close to the zero-temperature 
kinetic Ising model studies of ID chemical reactions. ( 8 > 13 ' 14 ) Detailed results are obtained 
in the appropriately defined continuum limit of the discrete dynamics (Section 5). These 
results are analyzed (Section 6) with emphasis on the length scales associated with the two- 
point correlation function. The structure-factor scaling at coexistence and its breakdown 
off coexistence are elucidated. 



2. DEFINITION OF THE MODEL 

In this section we define the model in ID. We also describe the extension to D > 1. 
However, the emphasis in this work is on the ID case, and the notation is introduced 
correspondingly. Thus, we consider spin variables aiit) = ±1. Time evolves in unit steps: 
t = 0,1,2,.... It is convenient to put the spin variables only at even lattice sites, i = 
0, ±2, ±4, ... for even times t = 0, 2, . . .. Similarly, for odd times t = 1, 3, ... we put spins 
at odd lattice sites i = ±1, ±3, — 

The values <7i(t + 1) for (t + 1) > will be determined stochastically by the dynamical 
rule incorporating interface diffusion and pairwise annihilation leading to cluster coarsening, 
and also spontaneous formation of the stable, +, phase from the unstable, — , phase thus 
attempting to model cluster growth in nucleation processes. If the "parent" spins, ai-i(t) 
and <7j_|_i(t) are both +1 or both —1, then the "offspring" is first set to +1 or —1, respectively. 
However, the —1 value is then overturned with probability p. If the "parent" spin values 
are opposite, the "offspring" is first set to one of them randomly. However, the —1 value is 
again overturned spontaneously with probability p. 

The first updating step (corresponding to setting p = 0) describes a symmetric "voter 
model" type dynamics. ( 1 ' 2 ' 6 ' 12 ) Indeed, the ordered, all + or all — , regions are unaltered. 
However, each interface between the neighboring + and — spin pairs hops one lattice spacing 
to the left or to the right with equal probability. On encounter, interfaces annihilate pairwise 
leading to cluster coarsening. Note that interfaces can be viewed as located at the odd sub- 
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lattice at even times and at the even sublattice at odd times. An important property of sym- 
metric voter models and related particle (here, interface) diffusion-with-annihilation models 
is the decoupling of the hierarchy of recursion relations for the correlation functions/ 1 ' 6 ) as 
well as the diffusive nature of the resulting equations in the continuum limit, which have 
allowed derivation of several exact results. ( 5_14 ) 

The second step of updating, i.e., the spontaneous spin-flips —1 — > +1, is introduced 
here as the means to break the ± symmetry and model formation of the + phase by growth 
from the — phase. Thus, we will be interested in the results for p <C 1. The key observation 
(not limited to D = 1) is that the correlation function hierarchy can be set up in such a 
way that decoupling properties reminiscent of the symmetric case are obtained. 

At each lattice site and for each time t > we introduce two random variables, Ci(t) 
that takes on values or 1 with equal probability, and Oi(t) which is with probability 
(1 — p) and 1 with probability p. The stochastic dynamics is defined by 



<Ti(t+l) = [1 - 6i(t + 1)] {d{t + lK_!(t) + [1 - (i{t + 1)] <T i+ l(t)} + 6i(t+l) (2.1) 

Given that all the random variables £ and are statistically independent, one can easily 
verify that the rule (2.1) correctly incorporates the dynamics as described in the preceding 
paragraphs. 

In calculating the averages, we can use the properties a 2 = 1, ( 2 = (, 9 2 = 9 at fixed 
time and lattice coordinate. Furthermore, Q(t) = §, @i(t) = p, where the overbars denote 
statistical averages. 

However, we still have to specify the initial values at t = 0. Either the values (Xj(0) 
can be given deterministically or quantities involving the later-time values Oi{t > 0) can be 
averaged over the distribution of the initial conditions. Here we prefer the latter option; 
we assume that the initial values are random and uncorrelated, with o"i(0) = ji. Thus, 
— 1 < /tx < 1 is the initial magnetization: each o"i(0) takes on values +1 and —1 with 
respective probabilities (1 + ji)/2 and (1 — /u)/2. 

Let us now define the average quantities that will be considered in this study. Firstly, 
due to translational invariance of the initial conditions and of the dynamical rule (2.1) after 
averaging over the random variables, the magnetization, m(t) = (Ti(t), depends only on time 
and satisfies the recursion (t > 0) 

m(t+ 1) = (1 -p)m(t) +p (2.2) 

with 

m(0) = fi (2.3) 
Similarly, the two-point correlation function, 



G n it) = ai(t)a i+n (t) (2.4) 

depends only on the distance n = 0, 2,4, . . . between the spins. The recursion relation for 
the two-point function for n > and t > is easily derived from (2.1): 
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G n {t + 1) = 1(1 - p) 2 [G n _ 2 (t) + 2G n (*) + G n+2 (t)] + 2p(l - p)m(t) + p 2 (2.5) 



where the initial and boundary conditions are 

G n=0 (t >0) = 1 and G n>0 (t = 0) = i? (2.6) 

The decoupling of the equations for the correlation functions should be obvious at this 
stage. Due to the linearity of (2.1) in a, the /c-point averages at (t+1) are determined by the 
k, k — 1, . . .-point averages at time t. This property is further amplified for the connected 
correlation function defined by 

C n (t) = G n (t)-m 2 (t) (2.7) 
Indeed, the appropriate recursion is 

C n (t + 1) = 1(1 - pf [C n _ 2 (t) + 2C n (t) + C n+2 (t)} (2.8) 

so that the m(t)-dependence enters only via the boundary conditions. The equivalents of 
relations (2.6) are 

C n=0 (t > 0) = 1 -m 2 (t) and C n>0 (t = 0) = (2.9) 

where we used (2.3). 

The probability to find an interface at i, in the interstice between the two spins <7i±i, 
is given by 

Pit) = 1 [1 - G 2 (t)} = 1 [1 - m 2 (t) - C 2 (t)} = 1 [C (t) - C 2 (t)} (2.10) 

Similarly to m(t), this quantity is translationally invariant. Both m and p can be also 
considered as the order-parameter and interface-number densities if we allow for the fact 
that they are defined per site of the lattice of twice the spacing of the original ID linear 
system of sites labeled by i. In fact, all our calculations will be with dimensionless variables 
such as distance n and time t. One can of course introduce dimensional length and time 
scales which has been a common practice especially in the continuum limit. However, 
we found that no new useful physical insight in gained, while the equations become more 
complicated. Thus, we use the dimensionless variables throughout. 

Before we outline the extension to D > 1, which will be detailed elsewhere, let us 
emphasize three appealing features of the ID model: the property that only two parent 
spins "vote" at each time step, the linearity of the evolution rule (2.1), and the fact that for 
p = there are already many results available, in particular, the relation to the interface 
motion and the interpretation of cluster coarsening due to interface annihilation. 

The simplicity of two-spin voting and the linearity of the dynamical rule (implying, 
essentially, solvability) can be extended to D > 1 by using the idea of updating along 
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different axes in each time step.^ 2 -* Consider spins <Ji 1 i 2 ...i D (t). For time steps t = — > 1, 
D — > D + 1, 2D — > 2D + 1, • • ., the rule (2.1) is used along axis 1, i.e., with ii varied as in 
(2.1) while ii-, - ■ -i^d kept the same on both sides of the relation. Similarly, for time steps 
1 — > 2, D + 1 — > D + 2, . . ., the update relation involves the index Z2 along axis 2, and so 
on. In D time steps, the cycle of the axis indices is complete. 

Regarding the availability of exact results for p = and the interpretation of the 
dynamics of the broken bonds connecting ± spin pairs, the D > 1 results, see Ref. 10 and 
literature cited therein, are understandably less numerous than those available for D = 1. 
In fact, it has been argued^ 1 ' 15-17 -* that symmetric voter-model type dynamics cannot lead 
to cluster coarsening in D = 3 and higher. In D = 2 the clusters do grow^ 1 ' 15-17 - 1 but 
the process can no longer be described by a simple cluster-size scaling. ^ 1,2 ' ) Quite generally, 
many open questions remain for D > 1. 



3. GENERATING FUNCTION FORMULATION 

The recursion relation (2.2) for the magnetization is trivial to solve, 

m(t) = l-(l-fi)(l-pY (3.1) 

However, the solution for the correlation function can be obtained in a simple form only in 
terms of the generating functions 

CO 

B n (v) = J2v t C n (t) (3.2) 



t=0 



Indeed, relations (2.8) yield, for n = 2,4, . . ., 



B n = ^(1 - p f (B n _ 2 + 2B n + B n+2 ) (3.3) 



while for n = we get 



B0{V) ~ l-(l-p)v - l-(l-pYv (3 ' 4) 

where we used the conditions (2.9) in deriving both (3.3) and (3.4). 

The second-order difference equation (3.3) has two linearly independent solutions of the 

form 

B n (v) oc b n/2 (3.5) 

where b{v) is a root of the quadratic characteristic equation. However, one can check that 
only one of the two roots yields the solution which converges exponentially to zero as n — > oo 
(for fixed v is the vicinity of 0). The other root yields exponentially divergent terms. After 
some algebra we arrive at the expression 
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B n (v) =B (v) 



1 - y/1 - (l~p) 2 V 

(1 — p) 2 v 



i/2 



(3.6) 



where n = 2, 4, . . .. 

The result (3.6), when expanded in powers of v, yields C n (t) as the t th Taylor series 
coefficient. However, the expressions thus obtained involve double sums and are rather 
unilluminating. The continuum limit results derived in Sections 5 and 6 provide a more 
useful source of physical insight on the nature of the dynamics. 

Our main interest presently will be in the expression of the generating function for the 
interface density p(t), see (2.10). This quantity is the t th Taylor coefficient of the function 



1 



B (v)-B 2 (v) 



(3.7) 



Explicit calculation yields the result 



pit) = 2(1 -p)(l- pY [1 -S t (l- p)] - (1 - pY(l - p)« [1 - S t (l)] 
where we defined the finite sum 



(3.8) 



St(a) = 



(2k)\a 



k=0 



kl(k+l)\2 2k + i 



Note that 



(3.9) 



Soo(a) = (l - Vl - a) I a 
where the t = oo Taylor series converges for all a in [0, 1] . 



(3.10) 



4. DIFFUSION AS OPPOSED TO SPIN-FLIP 

It is instructive to consider two models which represent the extremes of diffusion only 
or no diffusion at all, as far as interfacial dynamics is concerned. If we set p = in our 
model, the only processes are those of interface hopping and pairwise annihilation. Thus, 
the model falls in the class of the diffusion-limited chemical reactions A+A^inert, which, 
as well as related models, were studied extensively ( 5 ~ 14 ^ in ID. The density of interfaces 
reduces to 

t k , (k %^ <P = 0) (4.D 

=t+l V ' 



P(0) 



where generally, 



p(0) = (l-^ 2 )/2 
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The density of interfaces decreases monotonically and smoothly for discrete time steps 
t = 0^1^2^ — For large times, we have 

p(t) = 2p(0)/^rt (p = 0) (4.3) 

where the ~t -1 / 2 law is consistent with the previous exact calculations for these reactions. 

The other extreme would be to have no diffusion at all. For this, however, we have to 
modify our model. Thus, let us consider a model of ±1 spins with the only dynamical process 
consisting of spin-flips —1 — > +1 with probability p. Since the spins are uncorrelated, the 
model is trivial to solve. Indeed, the dynamical equation (2.1) is replaced by 

<n(t + 1) = [1 - 6i(t + 1)] <n(t) + 6i(t + 1) (4.4) 

where we now assume that the spins are located on the even sublattice at all times. 

The magnetization obeys the same equation (2.2), which simply reflects the fact that in 
our more complicated model diffusion of interfaces conserves the order parameter. However, 
in the new, uncorrelated-spin model (USM), all the /c-point correlations factorize trivially, 
and as a result the connected correlations vanish identically (for distinct k coordinates). 
Specifically, we get 

p USM (t) = (1 - m 2 )/2 = {l-n)(l- pY - i(l - ^) 2 (1 - p) 2t (4.5) 

C^ M (t) = 5 n , [l-m 2 (t)] (4.6) 

The time-dependent length scale of interest in cluster growth is the average size of the 
dominant, +, clusters. More generally, one may consider the cluster size distribution, which 
was not obtained analytically. (For some asymptotic results in the diffusion-only model see 
Ref. 6.) One measure of the cluster size is p _1 (t). For the diffusion-only model this cluster 
size measure grows according to ~ \ft for large times. For the USM, it grows as ~ (1 — p)~ l 
(assuming < p < 1). 

Note, however, that this quantity is related to the short-distance properties of the two- 
point correlations, see (2.10). The moment or decay-tail definitions of the "correlation" 
length scales (defined in Section 6) are typically used to probe the fixed-time large-distance 
behavior of the two-point correlations in strongly fluctuating systems. The various length 
scales are not necessarily related. For the diffusive model (symmetric, phase coexistence 
case, p = 0), it turns out that all the length scales behave according to n ~ \ft (Section 
6) . The USM example is instructive as the opposite extreme: the two-point correlations are 
zero-range; see (4.6). However, the cluster size measure p _1 is well defined and diverges as 
t — > oo. The length scale properties will be further explored in Section 6. 



5. CONTINUUM LIMIT 

The continuum limit has been the standard framework for writing phenomenological 
equations in cases which are not exactly solvable, or where the precise microscopic dynamics 
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is not known or specified. If fact, the continuum limiting description provides a useful guide 
for the identification of the "universality classes" or at least general classes of models with 
similar properties. A simple-minded continuum limiting procedure would amount to the 
assertion that for t 1 and n ^> 1 the discrete variation can be replaced by smooth 
functional dependence on t and n. Formally, one then uses the expression 



f(t + At, n + An) = \exp 



XAt 3i 



exp 



'a A d 

AAn— 

On 



f(t,n) 



(5.1) 



A=l, A=l 



to expand in the derivatives, which are presumably small. The order of the expansion is 
conveniently monitored by collecting powers of A and A before setting these variables to 1 
in the final expressions. 

If we apply this procedure to the equation (2.2) for m(t), and keep the leading t- 
derivative, we obtain the equation 



with the solution 



dm 

— — = —pm + p 

at 



m(t) = !-(!- n) e - pt 



(5.2) 



(5.3) 



where we used (2.3). However, this result differs from the exact expression (3.1), which is, 
in fact, perfectly well defined for all real t > 0. The source of the difficulty is clearly in that 
the t-derivatives are small only for p <^1. This is an illustration of the well-known property 
that the continuum approximation can be used only in a limited part of parameter space. 

A better controlled procedure is to use properly rescaled variables so that the parame- 
ters A and A in the equivalent of (5.1) are actually small. For our problem, we set 



t = pt and A = p 



(5.4) 



x = y/pn and A = ^fp (5-5) 

where the t-rescaling is suggested by our consideration of m(t), while the n-rescaling is 
implied by the diffusive combination n/t 2 which is expected to survive the p — > limit. 

In terms of the new variables, the relation (2.8) for the two-point function has the 
leading terms in order p. The "continuum limit" two-point function C{x,r) satisfies the 
relation obtained by collecting these terms, 



dC n „ d 2 C 
— = -2C+ — 
or ox 



-2C+^ (*>0) (5.6) 



where the initial and boundary conditions (2.9) are replaced by 

C(x = 0,t > 0) = 1 -m 2 (r) and C(x > 0, r = 0) = (5.7) 

with 
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m(r) = 1 — (1 — n)e T 
The interface density in the continuum limit is approximated as follows, 



dC(x,r) 
dx 



(5. 



(5.9) 



x=0 



The reader should keep in mind that the continuum limit is an approximation valid 
asymptotically for < p C 1, t 1, n 1. The results must be properly interpreted. 
For instance, if taken literally, relation (5.9) would imply that the interface density is infinite 
at t = because the initial conditions for C(x, r) are step-like. In fact, the divergence is in 
the regime where the continuum limiting approximation breaks down; see the next section. 
The rescaling (5.4)-(5.5) also obscures the p = case. Indeed, the results must be properly 
expressed in terms of the variable x/r 2 = n/t 2 before taking the limit p — > 0. If fact, p = 
is reminiscent of the "critical-point" limit in which there are no small parameters to rescale 
n and t. Instead, only their "scaling combination" enters in the continuum limit. 

The solution of the equation (5.6) with conditions (5.7) is obtained by the Laplace 
Transform method. We omit the mathematical details and only quote the final expression, 







2(1-//) (l-(i)'< 
u) + l uj + 2 



- x 



(5.10) 



which inverse-transforms to 



C(x,t) 
where 



(5. 



1) 



erfc(a) = -= / e~ p dp 

\/7T 



(5.12) 



is one of the standard error functions, the properties of which are well known. Thus, 
the expression (5.11) can be used to analyze various properties of connected two-point 
correlations. Some such results will be presented in the next section. 



6. LENGTH SCALES AND BREAKDOWN OF SCALING 

Let us consider the large-x behavior of C(x, r) for fixed r > (and /i ^ 1). All three 
terms in (5.11) then follow the asymptotic large-argument behavior of the error function. 
The result turns out to be 



C(x — > oo, r) oc e 



-2r 



X 



exp 



x 
4t 



(6.1) 
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where we omitted the proportionality constant. The decay-tail length scale, n ta ii, is thus 
determined by the dependence on the diffusional combination x 2 /r = n 2 /t, 



"-tail ~ Vt (6.2) 
Consider next the moment-definition length scales. We define the k th moment, 

OO 

M fc (r) = J x k C(x,r)dx (6.3) 
o 

and the associated time-dependent length, nf-(t), 

^n k = (M k /M ) 1/k (6.4) 

In the evaluation of M k , the contribution due to the first term in (5.11) can be always used 
in its large-argument form, while the third term is originally a function of the diffusional 
combination (times e~ 2r ). The second term, however, can be written in such a diffusional- 
scaled form only for 



x 



>aiv / r + a 2 r (6.5) 



where from now on the coefficient notation cij will be defined to stand for "a slow varying 
function of r, of order 1, possibly /c-dependent (when implied by context)." The diffusional 
contribution to the moments is 

Mf ff) = a 3 r^ +1 )/ 2 e" 2T (6.6) 

In the range of smaller x, not satisfying (6.5), the error function in the second term in (5.11) 
becomes of order 1. In fact, the fixed- ir, large-time behavior 

C(x,t -> oo) oc e~ T ~ x (6.7) 

is explicit in the Laplace-transformed form (5.10) due to the rightmost pole singularity at 
uj = — 1. The added contribution due to this exponential behavior is of the form 

a 1 ^/r+a 2 T 

M^ exp) = a 4 e" T J x k e~ x dx (6.8) 
o 

For small r, the intergation will yield the same power of r as in (6.6). Thus, the moment 
length scales behave according to 

n k ~ ( T (*+D/2/ T i/2) 1/k p -i/2 = y fi {t <<; 1/p) (6 9) 

However, as r increases, the integral in (6.8) saturates at a value of order 1. Since the 
remaining time dependence, e _r , dominates that of the diffusive contribution (6.6), the 
length scales saturate at 
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n k ~ 1/y/p (t > l/p) 



(6.10) 



The crossover between the limiting behaviors occurs at t ~ l/p and is difficult to evaluate 
in closed form. 

We next turn to the density of interfaces and the associated length scale. A direct 
calculation of the right-hand side of (5.9) yields 

-4= - ^£e- 2r + 2(l-^)e- T erf(v^) (6.11) 

x/p V^t 

where we kept the approximation sign to emphasize that this result applies only for t 3> 1 
(as well as p < 1). Note that erf (a) = 1 — erfc(a). The limit p = is thus correctly 
reproduced; see (4.3). For small r, the first term in (6.11) dominates, and the associated 
length scale behaves according to 

Tipocp -1 -^ (*<l/p) (6-12) 

However, for large r the second term takes over. Noting that the function erf (a) approaches 
1 for large a, we conclude that 



n p ~ e^/y/p (t > l/p) (6.13) 

In the theories of structure- factor scaling/ 4 -' where the structure factor is defined as 
the spatial Fourier transform of C(n,t), assuming continuous coordinate n and time t, the 
momentum, q, dependence is scaled in the form n(t)q. In the direct-space notation this 
amounts to assuming that the coordinate dependence enters via n/h(t). It is tempting 
to associate h(t) with a typical cluster size measure. In practice, n is determined as the 
inverse of some momentum scale found at low or fixed g-values, corresponding to large or 
intermediate coordinate values n. 

Our results support this picture at coexistence , i.e., at p = 0. Indeed, due to the 
critical-point-like scaling expressed by the diffusive scaling combination n 2 /t, all length 
scales defined at short or large distances are essentially the same. The identification n ~ \ft 
is unambiguous. However, explicit expressions obtained for p > indicate two difficulties 
with the structure-factor scaling when the growth of the stable phase occurs off coexistence. 
Firstly, the identification of a unique length scale is no longer possible for large times for 
which the cluster size distribution deviates significantly from the symmetric case. All three 
length scales estimated behave differently for large t. Secondly, the two-point correlation 
function no longer has simple scaling properties. In fact, a more general conclusion, alluded 
to in Section 4, is that in such cases the length scale n p (t) is the appropriate one to use as 
a typical + cluster size. However, it is characteristic only of the short-distance coordinate 
dependence of the two-point function. 

In summary, we presented a solvable ID model of cluster growth. Our results indicate 
that the ideas of structure-factor scaling apply only to cluster coarsening at coexistence. 
Off coexistence, a typical stable-phase cluster size measure reflects only the short-distance 
properties of the two-point correlations; the full correlation function no longer obeys scaling. 
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